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ABSTRACT 

Cosmological shock waves are a ubiquitous consequence of cosmic structure 
formation. They play a major role in heating baryonic gas in the large scale 
structure of the universe. In this contribution we study the shock-heated gas in 
connection with shocks themselves, using a set of N-body/hydrodynamic simula- 
tion data of a ACDM universe. The distributions of shock speed and temperature 
of shock-hated gas both should reflect the depth of gravitational potential of as- 
sociated nonlinear structures, so their morphology. From their distributions we 
find, in addition to hot gas in and around clusters/groups and warm- hot inter- 
galactic medium (WHIM) with T = 10^ — 10^ K mostly in filaments, there is a 
significant amount of low temperature WHIM with T < 10^ K distributed mostly 
as sheet-like structures. The low temperature WHIM was heated and collision- 
ally ionized mainly by shocks with Vsh ^ 150 km s~^, while photo- ionization by 
the UV and X-ray background radiations is important for metal ions. Not only 
the WHIM with T = 10^ - 10^ K but also the WHIM with T < 10^ K make 
up a significant fraction of gas mass, implying the low temperature WHIM could 
be important in resolving the missing baryon problem. The shock-heated gas in 
filaments and sheets are manifested best through emissions and absorptions in 
soft X-ray and far UV. We confirm that the WHIM with T = 10^ - 10^ K makes 
significant contributions to the soft X-ray background, absorptions of highly ion- 
ized species such as OVII and OVIII in AGN spectra, and line emissions from 
OVII and OVIII ions, as pointed by previous studies. However, the WHIM with 
T < 10^ K is the major contributor to absorptions of lower ionized species such 
as OV and OVI, because these photo-ionized ions are most abundant in sheets 
of low density and temperature. On the contrary, lines of OV (630 A) and OVI 
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(1032 A) are emitted mostly from the WHIM with 10^ < T < 10^ K, because 
they are from coUisionally excited ions. 

Subject headings: cosmology: theory - diffuse radiation - large-scale structure of 
universe - methods:numerical - quasars: absorption lines 

1. Introduction 

In the current paradigm of the cold dark matter cosmology, the large scale structure of 
the universe has formed through hierarchical clustering of matter. Deepening of gravitational 
potential wells caused baryonic matter to accrete onto nonlinear structures. It was postulated 
that the baryonic matter behaved like a gas and "collisionless shocks" formed in the tenuous, 
magnetized, cosmic plasma via collective electromagnetic viscosities (Kang et al. 1996; Quilis 
et al. 1998; Miniati et al. 2000; Gabici & Pasquale 2003; Ryu et al. 2003). In fact, the 
gravitational energy of the baryonic gas associated with structure formation was dissipated 
into heat primarily via such shocks. The dissipated energy, then, governs the thermal history 
and nature of the gas in the universe, and manifests the large scale structure as well as its 
dynamics through emitted radiation (Kang et al. 1994; Cen et al. 1995; Cen & Ostriker 
1999a; Dave et al. 2001; Valageas et al. 2002; Furlanetto & Loeb 2003). 

In a recent work, Ryu et al. (2003) studied the properties of "cosmological shock waves" 
and their roles on thermalization of gas and acceleration of nonthermal particles in the large 
scale structure, using cosmological N-body/hydrodynamic simulations. The cosmological 
shocks were classified into two types. "External shocks" form around outermost surfaces 
that encompass nonlinear structures, so they are by nature accretion shocks that decelerate 
the never-shocked intergalactic gas infalling toward sheets, filaments and knots. "Internal 
shocks" are produced within those nonlinear structures by infall of previously shocked gas to 
filaments and knots, during subclump mergers, as well as by chaotic flow motions. External 
shocks have high Mach numbers of up to M ~ 100 due to low temperature of the accreting 
gas. Internal shocks, on the other hand, have mainly low Mach numbers of M ~ a few, 
because the gas inside nonlinear structures has been previously heated by shocks (external 
and/or internal) and so has high temperature. However, internal shocks are more important 
in energetics, because of higher preshock density. As a result, thermalization of gas and 
acceleration of cosmic rays occurred mostly at internal shocks. 

The intergalactic gas was heated mostly by cosmological shock waves. The shock-heated 
gas with T = lO"'' — 10'' K was studied previously, utilizing cosmological hydrodynamic 
simulations (Cen & Ostriker 1999a; Dave et al. 2001). This component, which is called a 
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warm-hot intergalactic medium (WHIM) , was estimated to contain ~ 30% of baryon mass 
at the present epoch and resides mostly in filaments with a median overdensity of ~ 10 — 30. 
Such WHIM can count for a significant fraction of the missing baryons at low redshifts 
(Fukugita et al. 1998; Fukugita 2003). X-ray emission from the WHIM may contribute 
to a fraction of the soft X-ray background radiation (XBR) according to Phillips et al. 
(2001). While there were recent reports of observations of the WHIM around clusters (see, 
e.g., Finoguenov 2003; Zappacosta et al. 2004b), detection of such emission still remains 
a challenging task due to low surface brightness. On the other hand, the WHIM may be 
better detected through absorption systems and emission lines in soft X-ray and far UV. For 
instance, the possibility of probing the WHIM with OVI absorbers was explored by Cen et 
al. (2001) and Fang & Bryan (2001). Absorption systems due to other Oxygen and metal 
ions in the intervening WHIM were studied by Hellsten et al. (1998), Chen et al. (2003) and 
Viel et al. (2003). Detections of such absorbers in X-ray observations including Chandra 
were reported (see, e.g., Nicastro et al. 2003; McKernan et al. 2003, among many). Emission 
lines from metal ions in the WHIM were studied by Yoshikawa et al. (2003), Furlanetto et 
al. (2003), and Fang et al. (2003). A possible detection of such emission hnes from a filament 
around Coma was reported (Finoguenov 2003). 

In this paper, we study the shock-heated gas in the large scale structure of the universe 
and evaluate its observational manifestations with data from a set of cosmological hydrody- 
namic simulations where shock waves were identified (see Ryu et al. 2003). Specifically, the 
properties of shock-heated gas are analyzed in connection with those of cosmological shock 
waves. From the shock speed and gas temperature distributions, we see that 1) the shocks 
with Vsh ^ 700 km s~^ and the hot gas with T > 10^ K are distributed in and around clus- 
ters/groups, 2) the shocks with 150 < Vsh < 700 km s"^ and the WHIM with T = 10^-10^ K 
are mostly in filamentary structures, and 3) the shocks with Vg^ ^ 150 km s^"*^ and the WHIM 
with T < 10^ K arc mostly in sheet-like structures. Here wc refer the uncollpased, collision- 
ally ionized, shock-heated gas as the WHIM and extend the low temperature bound of the 
WHIM to well below lO^K. While the WHIM with T = 10^ - 10^ K has been studied ex- 
tensively in many previous works, the "low temperature WHIM" with T < lO^K has not 
received much attention and will be the main subject of this study. 

In the next section, simulation data are described along with numerical details. The 
properties of shock-heated gas are described in §3, and its observational manifestations are 
discussed in §4. Summary is followed in §5. 
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2. Simulation Data 

We used the data from a set of N-body/hydrodynamic simulations of a ACDM universe 
with radiative coohng. The simulations were performed using a PM+Eulerian hydrodynamic 
cosmology code that was specifically designed to capture shocks with a high accuracy. It is an 
updated version of the code described in Ryu et al. (1993). For instance, the code now adopts 
the MC (monotonized central difference) limiter, instead of the original minmod limiter to 
achieve sharper resolution of discontinuities (see, e.g., LeVeque 1997). The ACDM model 
adopted the currently favored values of the following cosmology parameters: VLbm = 0.043, 
floM = 0.227, Oa = 0.73 {Qbm + ^dm + ftA = l),h = i/o/(100 km/s/Mpc) = 0.7, and 
(Tg = 0.8. These values are consistent with those fitted with the recent WMAP data (see, 
e.g., Bennett et al. 2003; Spergel et al. 2003). A cubic region of size L = 100/i^^Mpc at 
present was simulated inside the computational box with 1024^, 512^ and 256'^ grid zones for 
gas and gravity and with 512^, 256^ and 128^ particles for dark matter, allowing a "fixed" 
spatial resolution of Al — 97.7h~^ kpc — 390.6/i~^ kpc. We note that the simulations of 
different resolution have the exactly same large scale structure. The only difference is the 
Nyquist frequency in the initial condition, apart from spatial and mass resolutions. 

We adopted the tabulated coohng rate that had been calculated by following the non- 
equilibrium coUisional ionization of gas with a given metallicity, cooling from 10^'^K to 10^ 
K under the isobaric condition (Sutherland & Dopita 1993). The radiative coohng rate 
based on equilibrium coUisional ionization is not a good approximation for T < 10^ K, since 
recombination lags cooling and so ionization fractions depend on thermal and ionization 
history of gas. The non-equilibrium cooling rate for the gas cooling under the isobaric 
condition, however, becomes a function of temperature only, if the initial temperature is 
high enough to ensure the initial ionization equilibrium {e.g., T > 10^ K) and if only 
two-body coUisional processes are included (see, e.g., Sutherland & Dopita 2002). In our 
simulations, gas was allowed to cool only down to T ~ 10^ K and coohng was set equal to 
zero for T < lO^K. For void gas, the minimum temperature was set as the temperature of 
cosmic microwave background radiation, i.e., Tj^in — Tcmb{z). For metallicity, we adopted 
the mean metallicity-density relation, 

(1) 

which had been suggested by Cen & Ostriker (1999b), using the data from a cosmological 
simulation with feedbacks into the IGM from stars. The simulations are the same ones 
previously used in Ryu et al. (2003) for the study of cosmological shock waves, except 
radiative cooling is included. 
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A few points on numerical details are noted. 1) Our simulations did not include star and 
galaxy formations and feedbacks from those stars, galaxies and AGNs, nor photo-ionization 
heating of the intergalactic gas. This is because our primary focus is the heating by cos- 
mological shock waves and the properties of shock-heated gas in filaments and sheets. Non- 
gravitational feedbacks arc expected to be most important in high density regions inside 
clusters and groups. Through previous studies by Cen & Ostriker (1999a) and Dave et al. 
(2001), it is known that the intergalactic gas was heated primarily by cosmological shock 
waves. Feedbacks may have contributed only ~ 20% or so of the thermal energy of the 
intergalactic gas, so they should have played only minor roles in the thermal evolution of 
the WHIM in filaments and sheets. On the other hand, photo-ionization heating should 
have affected mostly diffuse gas in void regions and Ly a clouds, heating it to T ~ 10^ K 
(see, e.g., Valageas et al. 2002). But it should not have affected much the properties of the 
WHIM, since most of the WHIM has T > 10^ K in our simulations. When we estimated the 
amount of shock-heated gas, we excluded the components that would have been heated pri- 
marily by photo-ionization (see the next section for details). However, the photo-ionization 
of metal ions, such as Oxygen, due to the background radiation was taken into account 
in post-processing data analyses, since such photo-ionization dominates over coUisional ion- 
ization for the WHIM in filaments and sheets. 2) Because of the fixed-grid-based nature, 
our data have a limitation in studies of high density regions such as cores of clusters and 
groups of galaxies, compared to those using SPH or AMR codes (Kang et al. 1994; Frenk 
et al. 1999). However, the same nature of fixed grid-spacing turns out to be an advantage 
in studying the low density gas in filaments and sheets. The typical thickness of filaments 
and sheets is ~ l/i~^Mpc or larger (see Fig. 2 below). In addition, the internal shocks inside 
those nonfinear structures have a mean separation of ~ l/i~^Mpc too (Ryu et al. 2003). As 
shown below, numerical resolution is important in studies of the gas inside those nonlinear 
structures. The fixed grid-spacing of AZ = 97.7/i~^kpc, ~ 1/10 of those scales, in our highest 
resolution simulation seems to be good enough to resolve filaments and sheets and shocks in 
those structures (Ryu et al. 2003). 3) In order to establish the connection between shock- 
heated gas and cosmological shocks, we identified shock surfaces by the procedure described 
in Ryu et al. (2003). Only the shocks with Vgh > 15 km s~^ were identified, although the 
code itself captures shocks of any speed. 
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3. 



Shock-Heated Gas in Filaments and Sheets 



3.1. Temperature and Density of Shock-Heated Gas 



We begin our discussion with the temperature and density of shock-heated gas in the 
simulated universe. Figure 1 shows the mass distribution in the plane of gas density and 
temperature z — from different simulations. The purple line along T ~ 10^ K in high 
density regime in the simulations with 1024^ and 512^ grid zones is an artifact of numerical 
simulation. In order to isolate the shock-heated WHIM, we draw two lines in the figure: 1) 
The straight line represents the "Lyman-a equation of state" of Valageas et al. (2002), 



which was adjusted for different values of cosmological parameters employed (see also Hui 
& Gnedin 1997, among many references). The line comes from the interplay of radiative 
heating by the UV background with the expansion of the universe. The region of overdensity 
5 < 10 that hes below this line would have been heated primarily by the background UV 
radiation to the temperature depicted by the line. 2) The curved line in the high density 
regime delineates the "fast-cooling" dense region, which was adopted also from Valageas et 
al. (2002) again with a suitable modification due to different values of parameters. It comes 
from the equality of the effective cooling time due to radiative cooling plus heating by the 
background UV radiation and the Hubble time. The high density region [S > 10) below the 
Lyman-a line and right to the fast-cooling line would have been collapsed into objects such 
as stars and galaxies or would be in the form of bound objects, if resolution were not limited 
by fixed grid-spacing in our simulations (see e.g., Furlanetto et al. 2003, for the distribution 
of gas from an SPH simulation). Since our focus is on the intergalactic gas that was mainly 
heated by shock waves, we exclude in our discussion and analyses below the photo-heated 
gas as well as the gas in collapsed/bound objects. 

Several points are apparent in Figure 1. 1) In addition to the hot component with 
T > 10^ K and the WHIM with T = 10^ — 10^ K, which were discussed in previous studies of 
shock-heated gas (Cen & Ostriker 1999a; Dave et al. 2001), there exists a significant amount 
of low temperature WHIM in the range of T < 10^ K. The mass fraction of the gas in each 
component is the followings: the hot component ~ 2.2%, the WHIM with T = 10^ — 10'' K 
~ 24% and the the WHIM with T < 10^ K ~ 13% in our highest resolution simulation data 
(see Table 1). In the hot component, the gas within 1 /i~^Mpc of clusters/groups hotter 
than IkeV was not included. Most of the low temperature WHIM has the temperature of 
10^ K < T < 10^ K, while the mass of the WHIM with T < 10^ K is ~ 3% of total gas 
mass. 2) Our estimate of mass fraction for the WHIM with T — 10^ — 10^ K is comparable 




(2) 
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to, although a httle smaller than, ~ 30% from simulations using the PTreeSPH and AMR 
codes (Dave et al. 2001). A part of the difference would be due to feedbacks from stars, 
galaxies, and AGNs, which were ignored in our simulations. 3) It is interesting to note 
that the distribution of our WHIM follows very closely the "equation of state for the warm 
IGM", analytically suggested by Valageas et al. (2002). However, their estimation of mass 
fraction for the "warm IGM", ~ 24%, is somewhat smaller than our estimation for all the 
WHIM with T < 10^ K. 4) Gas density reaches to higher values and coohng becomes even 
more efficient with the dependence at higher resolution. The resulting trend is that there 
is more gas with high density and low temperature in higher resolution simulations. As a 
result, the mass fraction of the WHIM continues to decrease as resolution increases in our 
simulations (see Table 1). 



3.2. Spatial Distributions of Shock Waves and Shock-Heated Gas 

In this subsection, we investigate the spatial distribution of shock-heated gas in con- 
junction with that of shock waves themselves. As mentioned in the introduction, two types 
of shocks can be distinguished. "External shocks" are accretion shocks around nonlinear 
structures, whose speed is basically the infall velocity of accreting flows. "Internal shocks" 
are products of infall, merging, and chaotic flow motions induced by gravitational interac- 
tions inside nonlinear structures. So the shock speed of both types should reflect the depth 
of gravitational potential wells of nonlinear structures, and hence their morphology. In our 
simulations, 1) high-speed shocks with Vsh ^ 700 km s~^ are found mostly around and inside 
knot-like structures of clusters/groups, 2) shocks with 150 km s~^ < Vgh ^ 700 km s~^ mostly 
around and inside filamentary structures, and 3) low-speed shocks with Vsu ^ 150 km s~^ 
mostly around and inside sheet-like structures. The left panels of Figure 2 show the locations 
of external shocks, divided into the three speed groups, in a two-dimensional slice ai z = 
from the 1024^ simulation. The shocks encompass and outline nonlinear structures. The 
structures with width of ~ 1/i^^Mpc in the top-left panel are mostly sheet-like structures 
with shallow potentials, while the thicker structures with width of > 2/i^^Mpc in the mid-left 
panel are filaments (see also Fig. 3 below). The structure in the bottom- left panel contains 
a cluster with X-ray emission-weighted temperature, ~ 2.4 keV, and X-ray luminosity, 
Lx ~ 4.1 X 10^^ /i erg s~^ (note this cluster has a high luminosity for its temperature due to 
excessive cooling without feedbacks from stars and galaxies). 

The morphological distinction can be revealed more clearly in three-dimensional volume 
renderings for the distribution of shocks waves. The image in Figure 3(a) displays all shocks 
(external and internal) inside the full box of (100 /i^^Mpc)^. At a quick glance, filamentary 
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structures look dominant. However, a careful inspection, especially near the edges of the 
box where structures do not overlap, exhibits lots of of sheet-like structures. In order to 
take a close look at sheet-like structures, we show in Figures 3(c) and (d) the distribution 
of low speed shocks with Vsh < 150 km s~^ in a region of (25 h~^Mpc)^, which is marked 
in the top-left panel of Figure 2. For filamentary structures we show in Figures 3(e) and 
(f) the distribution of shocks with 150 km s~^ < Vgh < 700 km s~^ in a region of (31 
h~^Mpc)^, which is marked in the mid-left panel of Figure 2. For further reference, volume 
rendering movies for shocks of different speed groups in the two regions shown in Figure 3(c) 
to (f) as well as in the region covering the complex in the bottom-left panel of Figure 2 are 
posted at "http://canopus.chungnam.ac.kr/ryu/shock.html". The figures (along with the 
movies) demonstrate that the morphology of nonlinear structures can be revealed through 
the distributions of shocks of the three speed groups. 

Naturally, it is expected that the spatial distribution of shock-heated gas is closely 
related to that of cosmological shock waves. Figure 3(b) shows a three-dimensional volume 
rendering image for the temperature of shock-heated gas in the same perspective as that of 
Figure 3(a). The images confirms that the two distributions are indeed very similar. The 
right panels of Figure 2 show the gas distribution in three temperature ranges, 10^ < T < 
10^ K, 10^ < T < 10^ K, and T > 10^ K, in the same two-dimensional slice as in the 
left panels. With the choice of temperature ranges, different morphology defined by shock 
waves is revealed once again by gas temperature too: i.e., the hot component mostly in 
knot-like structures, the WHIM with T — 10^ — 10^ K mostly in filaments, and the low 
temperature WHIM with T < 10^ K mostly in sheet-like structures. Hence our results 
suggest the existence of a component of the intergalactic gas, the low temperature WHIM 
with T < 10^ K, that has not been explored in details in previous studies. 



4. Shock-Heated Gas in Observations 

4.1. X-Ray Emission from Shock-Heated Gas 

As noted in the introduction, it was suggested that the X-ray emission from the WHIM 
may contribute significantly to the soft XBR. Phillips et al. (2001) estimated that the emis- 
sion would contribute ~ 10% of the observed XBR in the 0.5—2 keV range. They used ACDM 
simulation data in which radiative cooling, feedbacks from galaxies and photo-ionization were 
included, and 512^ grid zones were used for the box of 100/i~^Mpc size. On the other hand, 
based on analyses of observational data, Zappacosta et al. (2004a) suggested that most of 
the soft XBR could come from the WHIM. In order to asses how numerical details includ- 
ing finite resolution would affect the prediction, we calculated the X-ray emission from the 
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shock-heated gas outside clusters/groups using our simulation data with 1024^, 512^, and 
256^ grid zones. Similarly as Phillips et al. (2001), the gas in spherical regions of l/i~^Mpc 
radius around clusters/groups with T^. > 1 keV was excluded, since they would be identified 
as discrete X-ray sources and be removed from contributors to the background radiation. 
Also the gas below the Lyman-a and fast-cooling lines in Figure 1 was excluded. The X-ray 
spectrum code for an optically thin gas, MEKAL (Mewe et al. 1985), was used to calculate 
the mean proper volume emissivity, e{E,z), from the entire simulation box z — 2, 1.5, 
1, 0.5, 0.2, and 0. The mean background intensity at 2; = 0, J(-E'), was calculated then by 
integrating e(£', z) from 2; = 2 to 2; = 0. Also the energy flux, dF/d\ogT, in the X-ray band 
of 0.5 — 2 keV from the gas with temperature between logT and logT-l-iilogT was calculated 
from the proper volume emissivity, e{E,T,z), as in the calculation of J{E). The upper and 
lower panels of Figure 4 show J{E) and dF/dlogT, respectively, from the simulation data 
of different resolution. 

Two points are noticed in Figure 4. 1) From the bottom panel, we conflrm that among 
the shock-heated gas, it is the WHIM with T ~ 5 x 10^ — 10^ K that contributes most to the 
soft XBR. 2) The amount of X-ray emission increases with resolution, although the amount 
of shock-heated gas that emits such radiation decreases (see Table 1). It is because there 
is systematically more gas with higher density in higher resolution simulations. The X-ray 
emission in our two high resolution simulation (1024^ and 512^) differs by only a factor of 
two or so in £^ > 0.5 keV. The fractional contribution to the observed XBR (see PhiUips et 
al. 2001, for discussions on observations) is largest in the 0.5 — 1 keV range, and estimated to 
be ~ 30% in the 1024^ simulation and ~ 15% in the 512^ simulation. Our estimation from 
the 512^ simulation is comparable to, although a little larger than, ~ 10% of Phillips et al. 
(2001), which has the same spatial resolution. Although our simulations and that in Phillips 
et al. (2001) were performed with the same numerical code, different treatments of additional 
physical processes as well as details of numerical schemes should have led to the difference. 
However, our estimation from the highest resolution simulation with 1024'^ grid zones is ~ 3 
times larger than theirs, indicating that the WHIM could be a major contributor to the soft 
XBR. 



4.2. Absorption Systems of Shock-Heated Gas 

It was suggested that the properties of the shock-heated gas in filaments and sheets 
may be studied best by analyzing the soft X-ray and far UV "absorption systems" or "X-ray 
forests" in the spectra of distance quasars and AGNs. Especially, the OVI absorption systems 
in far UV and the OVII and OVIII absorption systems in soft X-ray were studied extensively 
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by several authors (see §1 for references). We also study the properties of such absorbers 
in our simulated local universe ai z = 0, following the procedure described below. First, 
we constructed a table for the fractions of Oxygen ions as a function of gas temperature 
and density. Since photo-ionization of Oxygen ions is important for the low density gas 
in filaments/sheets, the X-ray background radiation of Miyaji et al. (1998) and the UV 
background radiation of ShuU et al. (1999) were included, following Chen et al. (2003). The 
photo-ionization code, CLOUDY (Ferland et al. 1998), was used to compute the fractions of 
Oxygen ions in a grid of temperature and density. Figure 5 shows the resulting fractions of 
Oxygen ions at hydrogen number density —7 < lognn < — 4 as a function of temperature. 
For comparison, the fractions when no ionizing background radiation presents are also shown. 
It is obvious that the photo-ionization of Oxygen ions by the adopted background radiations 
is very significant in the temperature range that corresponds to filaments and sheets (T < 
10^ K). The most noticeable point is that, with the ionizing background radiations, the 
fractions of highly ionized Oxygen ions increase at lower density. According to Figure 1 the 
WHIM have the hydrogen number density, —6 < lognu ^ —5, so the fractions of OVI, 
OVII and OVIII ions are highest in those components. This is one of the reasons why these 
Oxygen ions were studied extensively in previous works. 

Next we calculated the column density distributions of absorption systems due to these 
Oxygen ions. In simulation data, clouds of shock-heated gas were identified along the Nip — 
3xNg {Ng the number of grid zones in one-direction) line paths of L = 100/i~^Mpc a,t z — 0. 
Again the gas in spherical regions of l/i~^Mpc radius around clusters/groups with > I keV 
as well as the gas below the Lyman-a and fast-coohng fines were excluded. The metalficity 
was assigned at each grid zone according to equation (1). The left panels of Figure 6 show 
the cumulative distributions of absorption system column density for OV-OVIII ions: 



where A'^abs is the number of absorption lines with column density, A'o ionsj greater than given 
values. As expected from Figure 5, the OVll and OVlll absorption systems are strongest. 
In the 1024'^ simulation data, dN/dz ~ 1 for both A^ovii and A^ovni > lO^^cm"^, which is in 
a good agreement with Chen et al. (2003) and Viel et al. (2003). The column densities of OV 
and OVI ions are smaller, so dN/dz ~ 1 for A^ov ^ 3 x lO^^cm-^ and A^ovi ^ 2 x lO^^'cm-^. 
We note the column density distributions in the 1024^ and 512^ simulation data are converged 
within a factor of two. 

The right panels of Figure 6 show the mass fractions of OV-OVIII ions due to the gas 
with temperature between logT and logT-|- dlogT, 
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which was normahzed with the total Oxygen mass Aio- We find that the WHIM with 

T = 10^ — 10^ K contributes to most of the absorption systems of OVII (~ 70%) and OVIII 
(~ 80%). Hence the observations of those absorption systems would explore filamentary 
structures, as pointed out in Viel et al. (2003). However, the low temperature WHIM with 
T < 10^ K is the major contributor to the OV and OVI absorption systems, accounting for 
~ 55% and ~ 60%, respectively. This emphasizes the importance of the low temperature 
WHIM in some observations, as noted by Bergeron et al. (2002) in the study of OVI systems 
in a quasar spectrum. 



4.3. Emission Lines from Shock-Heated Gas 

Observations of emission lines from the shock-hcatcd gas in filaments and sheets still 
remain technically challenging. However, it was suggested that detections of such emissions 
could be possible in future X-ray missions such as MBE (see http:/ /www.ssec.wisc.edu/baryons/) 
and DIOS (Ohashi et al. 2004). In this subsection, we study the Oxygen emission fines ra- 
diated by the WHIM in our simulation data. First, a table for the emissivity of OIII (698 
A), OIV (549 A), OV (630 A), OVI (1032 A), OVII (574 eV), and OVIII (653 eV) lines was 
constructed on a grid of gas temperature and density, using CLOUDY (Ferland et al. 1998). 
Collisional ionization as well as photo-ionization due to the UV and X-ray background ra- 
diations described in §4.2 were included. Figure 7 shows the resulting emissivity from the 
optically thin gas of hydrogen number density —7 < lognn < —4. For comparison, the 
emissivity when no ionizing background radiation presents is also shown. From this figure 
one can expect that for T > 10^ K the strongest fines from the gas with — 6 < log ^ ^5 
in filaments and sheets (see Fig. 1) would be the OVII and OVIII lines. However, for 
T < 10^ K, OV and OVI would produce emission fines stronger than OVII and OVIII lines. 

Using the emissivity table, the Oxygen line emissivity from the gas in simulation data 
was calculated at 2; = 0. As in the calculation of absorption column density, only the shock- 
heated gas was included and the same metalficity-density relation was adopted. The left 
panels of Figure 8 show the mean emissivity, jo unes, of OV - OVIII lines from the gas with 

temperature between logT and logT + dlogT. Overall the strength is greatest for the OV 
and OVI lines, followed by the OVII and OVIII lines. While virtually all emissions of the 
OVII and OVlll lines come from the higher temperature WHIM with 10^ < T < 10^ K 
(~ 95% for both lines), the OV and OVI lines are emitted mostly from the WHIM with 
10^ < T < 10^ K, ~ 75% and ~ 80%, respectively. Contrary to the OV and OVI absorptions 
which get significant contributions from the WHIM with T < 10^ K, the line emissions of 
OV (540 A) and OVI (1032 A) come mostly from the WHIM with T > 10^ K. Our estimates 
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of line emissions are converged again within a factor of two in the 1024^ and 512^ simulation 
data. 

Our results indicate that the exploration of shock-heated gas through emissions in soft 
X-ray, such as the MBE and DIOS missions, would pick up a fraction of the WHIM with 
10^ < T < 10^ K (~ 10% of total gas mass in our estimate), missing the lower temperature 
WHIM with T < 10^ K. A search mission using the OV and OVI line emissions in far UV 
would detect a larger fraction of the low temperature WHIM with 10^ ^ 10® K (~ 14% 
of total gas mass), complementing the proposed missions in the soft X-ray band. 

We also calculated the specific intensity of Oxygen lines, Iq unes, by integrating the 
emissivity /47r along columns of 100^~^Mpc at ^ = 0. Note that the Hubble expansion 
speed over the scale of 100/i~^Mpc is (l/30)c, and the angular size of a single grid of Al 
at a distance of L = lOO/i^-'^Mpc is Al/L radian. So Jq imcs corresponds effectively to the 
intensity measured over a band of width AA/ A ~ 1/30 with a beam of solid angle 3.4' x 3.4' or 
smaller, when computed with the 1024^ simulation data. We integrated along Nip = 3 x Ng 
different line paths in the simulation box. The right panels of Figure 8 show the cumulative 
fraction of line paths with Iq unes larger than given values. They show that the fraction of 
randomly-chosen hue paths with Iq lines > 10~^ erg cm~^ sr~^ would be ~ 10~^ for the 
OV and OVI lines, while it would be 2 — 3 x 10~^ for the OVII and OVIII lines. However, 
the fractions of Iq imcs > 10^^ erg cm~^ s~^ sr~^ are /a? ~ 1 — 2 x 10^^ for the OVII and 
OVIII lines and larger than those for the OV and OVI lines. This reflects the fact that while 
the OV and OVI lines arc stronger overall, the OVII and OVIII lines arc from the WHIM 
with T > 10® K and density higher than the mean density of shock-heated gas. 



5. Summciry 

Hierarchical clustering induces cosmological shock waves in the course of large scale 
structure formation in the universe (Ryu et al. 2003). The intergalactic gas was heated 
mostly by such shocks (Cen & Ostriker 1999a; Dave et al. 2001; Valageas et al. 2002). 
In this paper, we studied the properties, spatial distribution, and possible observational 
manifestations of the shock-heated gas in filaments and sheets with a set of simulation data 
of a ACDM universe from a grid-based N-body/hydrodynamic code (Ryu et al. 1993). The 
nature of fixed grid-spacing makes such data suitable to the study of the shock-heated gas 
in fllamcnts and sheets. Shock-heating was focused, so feedbacks from stars, galaxies and 
AGNs, and photo-ionization heating were ignored, although the photo-ionization of metal 
ions by the background radiation were taken into account in post-processing data analyses. 
We expect that exclusion of such processes should not have weakened significantly the main 
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results of this work, since the gas in filaments and sheets was heated mostly by cosmological 
shocks. 

The speed of cosmological shock waves reflects the depth of gravitational potential of 
the associated nonlinear structures, so their morphology. We saw that 1) the shocks with 
Vsh ^ 700 km s~^ are distributed mostly around and inside clusters/groups, 2) the shocks 
with 150 km < Vsh ^ 700 km are mostly around and inside filamentary structures, 
and 3) the shocks with Vgh ^ 150 km s~^ are mostly around and inside sheet-like structures. 
The distribution of the shock-heated gas should be closely related to that of shock waves. We 
found that the WHIM with T = 10^ — 10^ K is distributed mostly in filamentary structures, 
while the low temperature WHIM with T < lO'^K is mostly in sheet-like structures. The hot 
gas with T > lO^K resides on knot-like structures as the intra-cluster/intra-group medium 
and the low density medium around them. The amount of shock-heated gas was estimated 
as follows: the hot component with T > 10^ K (excluding the gas inside clusters/groups) 
accounts for - 2.2% of total gas mass, the WHIM with T = 10^ - 10^ K for ~ 24%, and the 
WHIM with T < 10^ K for ~ 13% in our highest resolution simulation. Thus our results 
indicate the existence of the WHIM with T < 10^ K, which was heated by shocks of low 
speed with Vgh ^ 150 km s~^ and is distributed as sheet-like structures. We suggest that the 
low temperature WHIM with T < 10^ K contributes to a significant fraction of the missing 
baryons at low redshifts (Fukugita et al. 1998; Fukugita 2003). 

It has been suggested in previous studies (sec §1 for references) that the shock-heated 
gas in filaments and sheets can be manifested through emissions and absorptions in soft 
X-ray and far UV. 1) We found that the soft X-ray emission from the shock-heated gas in 
our highest resolution data contributes to ~ 30% of the observed XBR in the 0.5 — 1 keV 
range. This is somewhat larger than that previously suggested by Phillips et al. (2001). 2) 
Column densities of Oxygen ions along randomly selected line paths in the computational 
box were calculated. Such column densities can produce the absorption systems or X-ray 
forests in the spectra of distant quasars and AGNs (see, e.g., HcUstcn et al. 1998). The 
photo-ionization of Oxygen ions by the background UV and X-ray radiations is important 
for the gas in filaments and sheets due to low density and temperature. As a result, the 
fractions of Oxygen ions depend strongly on the local gas density. We found that, for OVII 
and OVIII, the absorption systems associated with the WHIM with T — 10^ — 10^ K account 
for 70 — 80% of the identified absorption systems, while the low temperature WHIM with 
T < 10^ K contributes ~ 55 - 60% of the OVI and OV absorption systems. 3) Finally the 
emission fines from Oxygen ions were calculated. They are emitted mostly from coUisionally 
excited ions. We estimated that ~ 95% of the OVII line (574 eV) and OVIII line (653 eV) 
emissions come from the WHIM with 10^ < T < 10^ K, while ~ 75% of the OVI line (1032 
A) and 80% of the OV fine (630 A) are emitted by the WHIM with 10^ < T < 10^ K. 
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Hence, we conclude that absorption systems in UV due to species such as OV and OVI 
provide the best chance to detect the low temperature WHIM with T < 10^ K. 
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Table 1. Mass fractions of shock- heated gas 



Component 1024^ 512^ 256^ 

Hot with T > W K'^ 0.022 0.024 0.029 

WHIM with T = 10^ - 10^ K 0.24 0.37 0.43 

WHIM with T < 10^ K 0.13 0.13 0.10 



The gas within 1 h ^Mpc of clusters/groups hotter than IkeV was not counted. 
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Fig. 1. — Mass distribution of shock-heated gas in the plane of gas density and temper- 
ature at 2; = in simulations with different resolution. For comparison, the distribution 
from the 1024"^ simulation without cooling is also shown. Purple indicates the highest mass 
concentration, and the color scale was set arbitrary to highlight the distribution. The black 
straight line presents the "Lyman-a equation of state" and the black curved line isolates the 
"fast-cooling" region of Valageas et al. (2002) (see text for details). 



Fig. 2. — Left panels: Locations of external shocks at 2; = in a two-dimensional slice of (100 
/i~^Mpc)^ from the 1024^ simulation. The shocks are classified into three groups according to 
shock speed, Vsh < 150 km s~^, 150 < Vsh < 700 km s~^, and Vsh. > 700 km s~^, respectively, 
from top to bottom. Right panels: Spatial distribution of shock-heated gas at z = in the 
same slice in three ranges of temperature: the WHIM with T < 10^ K, the WHIM with 
T = 10^ — 10^ K, and the hot component with T > 10^ K, shown in the top, middle and 
bottom panels, respectively. In each panel purple represents to the gas with temperature 
close to the upper bound and green close to the lower bound, while blue in the middle. 



Fig. 3. — (a) Distribution of shock waves (external and internal) at 2; = in the full box of 
(100 h^^Mpc)^. Color represents logarithmically-scaled shock speed from Vsh = 15 km s~^ 
(yellow) to 1, 500 km s~^ (purple) and higher, (b) Temperature distribution of shock-heated 
gas at 2; = in the same box. Color covers logarithmically-scaled temperature from T — 
10^ K (yellow) to 10^ K (purple) and higher, (c, d) Distribution of shock waves (external 
and internal) with Vsh = 15 km s~^ (yellow) to 150 km s~^ (blue-green) at 2; = in a region 
of (25 /i~^Mpc)"^, marked at the top-left panel of Figure 2. Tow different perspectives are 
shown, (e, f) Distribution of shock waves (external and internal) with Vsh = 150 km s~^ 
(green) to 700 km s^^ (red) at 2; = in a region of (31 h^^Mpc)^, marked at the mid-left 
panel of Figure 2. Tow different perspectives are shown. All images are from the 1024^ 
simulation. 
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Fig. 4. — Upper panel: Mean background intensity in soft X-ray from shock-heated gas, 
integrated from 2; = 2 to 2; = 0, in simulations with different resolution. Lower panel: X-ray 
flux in the 0.5 — 2 keV energy band from shock-heated g 3jS clS db function of temperature, 
integrated from 2; = 2 to 2; = 0, in the same simulations. The gas in spheres of a radius 
1 h~^Mpc around clusters/groups hotter than 1 keV was excluded in calculating the X-ray 
radiation. 
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Fig. 5. — Fractions of OIII to OVIII ions as a function of temperature for the gas of hydrogen 
number density log(?T,j^) = —7 to —4. Colhsional ionization as well as photo-ionization by the 
UV and X-ray background radiations (Miyaji et al. 1998; Shull et al. 1999) were included. 
For comparison, the fractions of the ions in colhsional ionization equilibrium are also shown. 
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Fig. 6. — Left panels: Number of absorption lines per unit redshift by intervening gas clouds 
at 2; = with the column density of Oxygen ions larger than A^^o ions- Right panels: Fractions 
of OV - OVIII ions relative to the total Oxygen abundance in shock-heated gas at 2; = as 
a function of temperature. 
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Fig. 7. — Emissivity of OIII to OVIII ion lines as a function of temperature for the gas of 
hydrogen number density log(nH) = -7 to -4: OIII (698 A), OIV (549 A), OV (630 A), OVI 
(1032 A), OVII (574 eV), OVIII (653 eV). ]i/n]j is in units of lO'^^g cm^ s'^ Collisional 
ionization as well as photo-ionization by the UV and X-ray background radiations (Miyaji 
et al. 1998; ShuU et al. 1999) were considered. For metallicity, Z = O.IZq was set. For 
comparison, the emissivity of the ion lines in collisional ionization equilibrium is also shown. 
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Fig. 8. — Left panels: Mean emissivity of OV to OVIII lines from shock-heated gas at z = as 
a function of temperature. Right panels: Cumulative fraction of line paths with the specific 
intensity of Oxygen lines larger than Iq lines- The intensity was calculated by integrating the 
emissivity divided by An along a length of 100/;.~^Mpc at 2; = 0. 
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